The Journal of Biological Chemistry 

© 2005 by The American Society for Biochemistry and Molecular Biology, Inc. 


Vol. 280, No. 8, Issue of February 25, pp. 6933-6941, 2005 

Printed in U.S.A. 


Nidovirus Sialate-O-Acetylesterases 

EVOLUTION and substrate specificity of coronaviral and toroviral 
RECEPTOR-DESTROYING ENZYMES*[S] 


Received for publication, August 23, 2004, and in revised form, October 26, 2004 
Published, JBC Papers in Press, October 26, 2004, DOI 10.1074/jbc.M409683200 

Saskia L. Smits$, Gerrit J. Gerwig§, Arno L. W. van Vliett, Arjen LissenbergtH, Peter Brizall, 
Johannis P. Kamerling§**, Reinhard VlasakttWl, and Raoul J. de Grooti III 

From the $Virology Division, Department of Infectious Diseases and Immunology, Faculty of Veterinary Medicine, Utrecht 
University, 3584 CL Utrecht, and §Bijvoet Center for Biomolecular Research, Department of Bio-Organic Chemistry, 
Section of Glycoscience and Biocatalysis, Utrecht University, 3584 CH Utrecht, The Netherlands, the \\Department of 
Molecular Biology and %%Applied Biotechnology, Department of Cell Biology, University of Salzburg, 

A-5020 Salzburg, Austria 


Many viruses achieve reversible attachment to sialic 
acid (Sia) by encoding envelope glycoproteins with re¬ 
ceptor-binding and receptor-destroying activities. Toro- 
viruses and group 2 coronaviruses bind to O-acetylated 
Sias, presumably via their spike proteins (S), whereas 
other glycoproteins, the hemagglutinin-esterases (HE), 
destroy Sia receptors by de-O-acetylation. Here, we 
present a comprehensive study of these enzymes. 
Sialate-9-O-acetylesterases specific for 5-Af-acetyl-9-0- 
acetylneuraminic acid, described for bovine and human 
coronaviruses, also occur in equine coronaviruses and 
in porcine toroviruses. Bovine toroviruses, however, ex¬ 
press novel sialate-9-O-acetylesterases, which prefer the 
di-O-acetylated substrate 5-N-acetyl-7(8),9-di-0-acetyl- 
neuraminic acid. Whereas most rodent coronaviruses 
express sialate-4-O-acetylesterases, the HE of murine 
coronavirus DVIM cleaves 9-O-acetylated Sias. Under 
the premise that HE specificity reflects receptor usage, 
we propose that two types of Sias serve as initial attach¬ 
ment factors for coronaviruses in mice. There are strik¬ 
ing parallels between orthomyxo- and nidovirus biology. 
Reminiscent of antigenic shifts in orthomyxoviruses, ro¬ 
dent coronaviruses exchanged S and HE sequences 
through recombination to extents not appreciated be¬ 
fore. As for orthomyxovirus reassortants, the fitness of 
nidovirus recombinant offspring probably depends both 
on antigenic properties and on compatibility of recep¬ 
tor-binding and receptor-destroying activities. 


The first and arguably most critical step of any viral infection 
is the attachment of the virion to the host cell. The molecules to 
which virions bind comprise a diverse collection of cell surface 
proteins, lipids, and carbohydrates (1). Quite a number of en¬ 
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teric and respiratory viruses use sialic acid (Sia) 1 either as 
their sole receptor or as an initial attachment factor (2). 

Sialic acids, a generic name for a large family of 9-carbon 
negatively charged monosaccharides, typically occur as termi¬ 
nal residues of glycoconjugates. Of the variety of Sias that 
occur in nature, many result from differential O-acylation at 
C4, C7, C8, and C9 of the parental molecule 5-A r -acetylneura¬ 
minic acid (Neu5Ac) (2, 3). Sias make ideal receptors because of 
their accessibility and abundance, but their wide distribution 
also has a downside. High avidity binding of virions to nonper- 
missive cells, to soluble sialylated compounds, or to sialylated 
mucins of the mucus barrier all might lead to loss of infectivity. 
In addition, loss of infectivity might result from virus self¬ 
aggregation. To ward against such hazards, some enveloped 
RNA viruses have developed an elegant strategy, which allows 
high affinity yet reversible binding to Sia receptors; they spec¬ 
ify for virion-associated receptor-destroying enzymes (RDEs). 
So far, two types of viral RDEs have been identified: sialidases 
(neuraminidases) and sialate-O-acetylesterases. 

RDEs with sialidase activity are found in para- and ortho¬ 
myxoviruses (4, 5). For example, the orthomyxoviruses, influ¬ 
enza virus A and B, express two envelope glycoproteins: the 
hemagglutinin (HA), which mediates binding to Neu5Ac and 
subsequent fusion of viral and host cell membranes, and the 
neuraminidase (NA), which catalyzes the removal of Neu5Ac 
from sialoglycoconjugates. A functional balance between the 
antagonistic HA and NA activities is critical for efficient 
entry (5). 

The other type of RDE, sialate-O-acetylesterase, was first 
discovered in influenza C virus (6, 7). This virus expresses only 
one envelope glycoprotein species, the hemagglutinin-esterase 
fusion protein (HEF), in which all three functions (receptor 
binding, RDE, and membrane fusion activity) are combined (7, 


1 The abbreviations used are: Sia, sialic acid; aNeu4,5,9Ac 3 2 Me, 5-IV- 
acetyl-4,9-di-0-aeetylneuraminic acid a-methyl glycoside; BCoV, bovine 
coronavirus; BSM, bovine submaxillary mucin; BToV, bovine torovirus; 
D bb , distribution of synonymous substitutions; E, esterase domain; ECoV, 
equine coronavirus; F, fusion domain; HA, hemagglutinin; HCoV, human 
coronavirus; HE, hemagglutinin-esterase; HEF, hemagglutinin-esterase 
fusion protein; K B , synonymous substitutions per synonymous site; MHV, 
mouse hepatitis virus; NA, neuraminidase; Neu5Ac, 5 -N- acetylneura¬ 
mi nic acid; Neu5,9Ac 2 , 5-A r -acetyl-9-0-acetylneuraminic acid; Neu4,5Ac 2 , 
5-A r -acetyl-4-0-acetylneuraminic acid; Neu5,8Ac 2 , 5-lV-acetyl-8-0-acetyl- 
neuraminic acid; Neu5,7Ac 2 , 5-lV-acetyl-7-0-acetyl neuraminic acid; 
Neu5,7(8),9Ac 3 , 5-lV-acetyl-7(8),9-di-0-acetylneuraminic acid; pNPA, p- 
nitrophenyl acetate; PToV, porcine torovirus; R, receptor-binding domain; 
RCoV-SDAV, rat sialodacryoadenitis coronavirus; RDE, receptor-destroy¬ 
ing enzyme; S, spike; SARS, severe acute respiratory syndrome; HPLC, 
high pressure liquid chromatography. 
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8). The structure of the HEF spike, a homotrimer of a class I 
membrane glycoprotein, resembles that of HA of influenza 
virus A (9). Like the HA protein, the HEF protein is post- 
translationally cleaved. The N-terminal subunit, HEF1, forms 
the globular top part of the spike and contains the receptor¬ 
binding (R) and RDE esterase (E) domains; the C-terminal 
subunit HEF2, together with segments of HEF1, contains the 
fusion domain (F) and forms the large, membrane-anchored 
stalk (9). HEF specifically recognizes 5-N-acetyl-9-0-acetyl- 
neuraminic acid (Neu5,9Ac 2 ) receptors, via a binding site at the 
tip of the globular head region. The catalytic pocket of the RDE 
constitutes a second Sia-binding site. The RDE, a serine hydro¬ 
lase with a Ser-His-Asp catalytic triad, destroys Neu5,9Ac 2 
receptors by removing the 9-O-acetyl group (6, 10). 

An RDE with sialate-4-O-acetylesterase activity was recently 
discovered in infectious salmon anemia virus (genus Isavirus, 
family Orthomyxoviridae), a segmented negative-stranded RNA 
virus of teleosts (11). There is no obvious sequence identity be¬ 
tween influenza C virus HEF and infectious salmon anemia virus 
hemagglutinin-esterase protein (HE), except for the regions that 
contain the catalytic residues (12). However, influenza C virus 
HEF1 does share considerable identity with the nidovirus HEs 
(i.e. class I envelope glycoproteins expressed by toroviruses and 
certain coronaviruses) (13-16). 

Coronaviruses and toroviruses (separate genera in the fam¬ 
ily Coronaviridae, order Nidovirales ) are evolutionarily re¬ 
lated, enveloped ( + )-strand RNA viruses of mammals and birds 
(17-19). All toroviruses identified so far carry an HE gene (20), 
but among coronaviruses, only those belonging to (sub)group 2 
do. The latter are exemplified by mouse hepatitis virus (MHV) 
and bovine coronavirus (BCoV). In SARS-associated human 
coronavirus (HCoV-SARS), which is an early phylogenetic 
split-off of group 2 (21), and in group 1 and 3 coronaviruses, 
represented by porcine transmissible gastroenteritis virus and 
avian infectious bronchitis virus, respectively, the HE gene 
is absent. 

The role of the nidoviral HE proteins and their importance 
for viral replication is still subject to debate. In the case of 
BCoV, infection in the presence of esterase inhibitors results in 
reduced virus titers (22), whereas HE-specific antibodies neu¬ 
tralize the virus (23). Other observations, however, indicate 
that the HEs are nonessential for replication in tissue culture 
cells (24, 25). In fact, during propagation in vitro of murine 
coronaviruses, HE expression is readily lost (13, 26). 2 More¬ 
over, infection experiments with such HE-deficient MHV labo¬ 
ratory strains would suggest that, at least for these coronavi¬ 
ruses, the HE protein is also dispensable in vivo. 

Entry of coronaviruses and toroviruses [i.e. attachment and 
subsequent fusion of viral and host cell membranes) is medi¬ 
ated by the spike envelope glycoprotein S (27, 28). MHV, HCoV- 
SARS, and group 1 coronaviruses, transmissible gastroenteri¬ 
tis virus, and HCoV-229E use specific cell surface glycoproteins 
as their main receptor (29-32), as is probably the case for all 
other coronaviruses and toroviruses. However, there is accu¬ 
mulating evidence that coronaviruses in addition use Sia re¬ 
ceptors (33-39). Contrary to what might be expected from their 
similarity to influenza C virus HEF or from their name, the 
coronaviral HEs seem to play only a minor role in viral attach¬ 
ment to Sias. Rather, the spike proteins appear to be the major 
Sia receptor-binding proteins (35, 38, 39). Equine torovirus 
Berne, which lacks a functional HE gene, induces hemaggluti¬ 
nation (40), suggesting that toroviruses utilize Sia as an at¬ 
tachment factor as well and again implicating the S protein as 
the main candidate for Sia receptor binding. 


2 A. Lissenberg and R. J. de Groot, unpublished observations. 


For all group 2 coronaviruses studied so far, the Sia receptor 
binding specificity of S matches the substrate preference of HE. 
For instance, BCoV and HCoV-OC43 specifically bind via their 
S proteins to glycosidically bound Neu5,9Ac 2 moieties (35, 39) 
and, in accordance, encode HEs with sialate-9-O-acetylesterase 
activity (22, 41). Likewise, MHV strain S binds to 5-N-acetyl- 
4-O-acetylneuraminic acid (Neu4,5Ac 2 ) receptors apparently 
via its S protein and encodes an HE with sialate-4-O-acetyles- 
terase activity (38). Sialate-4-O-acetylesterases are also pres¬ 
ent in the closely related rodent coronavirus strains MHV-JHM 
and rat sialodacryoadenitis coronavirus (RCoV-SDAV) (16, 38, 
42, 43). 

Toroviruses and the group 2 coronaviruses apparently ac¬ 
quired their HE genes from an unknown source via independ¬ 
ent heterologous RNA recombination events (13-15). Conceiv¬ 
ably, the acquisition of an HE protein increased viral fitness by 
allowing optimal usage of specific Sias as attachment factors. 
Still, our understanding of the evolution, function, and sub¬ 
strate specificities of the nidovirus HEs is incomplete. Both the 
toro- and the coronaviral HE genes have diverged considerably 
(20, 38). Coronavirus HEs have been described, which in their 
sequence are clearly distinct from those analyzed so far and for 
which the substrate specificity remains to be determined. With 
respect to the specificities of the toroviral HEs, no information 
is available at all. Here, we present a comprehensive study of 
nidovirus HE proteins based upon a combination of phyloge¬ 
netic, biochemical, and molecular virological techniques. In 
addition to the sialate-9-O- and sialate-4-O-acetylesterases de¬ 
scribed so far, we report the identification of a novel type of HE, 
which occurs in bovine toroviruses (BToVs) and displays a 
preference for Sias with a di-O-acetylated glycerol side chain. 
We already noted that exchange of HE genes via homologous 
RNA recombination is common among toroviruses (20). We now 
show that naturally occurring rodent coronaviruses also ex¬ 
change their HE as well as their S genes to an extent not 
appreciated before. At least on one occasion, this must have 
resulted not only in an antigenic shift but also in a shift in HE 
substrate specificity. The recombinational exchanges of surface 
glycoprotein genes in nidoviruses are reminiscent of the reas¬ 
sortment of HA and NA genome segments in influenza A and B 
viruses. Parallels between orthomyxo- and nidoviral biology 
with respect to attachment, envelope glycoprotein function, 
and consequences thereof are discussed. 

MATERIALS AND METHODS 

Phylogenetic Analyses —HE sequences were extracted from the NCBI 
data base for BToV strains Breda (BRV), B145 and B150 (accession 
numbers Y10866, AJ575379, and AJ575380), porcine torovirus (PToV) 
strains P-MAR, P10, and P4 (AJ575363, AJ575366, and AJ575364), 
MHV strains A59, JHM, and 2 (NC_001846, AAA46442, and 
AF201929), PuCoV (AJ005960), RCoV-SDAV (AF207551), equine coro¬ 
navirus (ECoV)-NC99 (AY316300), BCoV-Mebus (U00735), HCoV- 
OC43 (OC43; M76373), PHEV (AY078417), and influenza C virus 
strains Miyagi/4/93 (C/MIY/4/93; BAB84727) and Johannesburg/66 (C/ 
JHG/66; S07412). In MHV-A59 and MHV-2, the reading frames of the 
HE genes are interrupted by termination codons and, in the case of 
MHV-2, also interrupted by a 4-nucleotide deletion. To allow the inclu¬ 
sion of these strains in phylogenetic analyses, their HE genes were 
reconstructed in silico according to the MHV HE consensus sequence: 
MHV-A59, TGA (nucleotides 43—45) > GGA; MHV-2, TAA (nucleotides 
292—294) —> CAA; and insertion of the sequence AAAC at nucleotide 
position 525. 

Discrepancies were noted between sequences deposited for MHV- 
DVIM in the NCBI data base, indicating that virus stocks have been 
switched. The DVIM HE nucleotide sequence of Sugiyama et al. 
(AB008939) is 100% identical to the HE gene of MHV-S, as determined 
in our laboratory, 3 but only 73% identical to the DVIM HE sequence of 


3 A. Lissenberg and R. J. de Groot, unpublished observations; Gen- 
Bank™ number AY771997. 
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Table I 

Oligonucleotide primers 


Oligonucleotide 

Nucleotide sequence (5'-3') 

Polarity 

Positions 0 

Template 

995 

ACTTATTATTTTGTTGAAATG 

+ 

-18-3* 

MHV-S 

996 

ACCAATATACCCTAAACAAGA 

- 

23959-23979“ 

MHV-A59 

1218 

CCACCATGTTGAGGATGAGGGTTCGTCCACC 

+ 

1-26* 

P-MAR 

1240 

CTAATAACTACTTAAACAAAAACT 

- 

1261-1284 6 

P-MAR, P10, P4 

1345 

GAGACACTATCTTTAG 

+ 

27284-27299“ 

P10, P4 

1383 

TAGCATTTGGATTAAGCATAGA 

- 

27319-27340“ 

B145, B150 

1394 

CGAATTGAACCCATAAATAACACCAGTGTC 

- 

681-700 6 

BCoV 

1433 

AAGTTTGAGTAGCCACTTATC 

+ 

26773-26793“ 

B145, B150 

1454 

TGAGGGTGTTAACTTTACGCCTTA 

+ 

301-324 6 

BCoV 

1455 

GTAAAGTTAACACCCTCATAAAAAAT 

- 

291-318 6 

BCoV 

1502 

TCGTTCAGATATCTCATCTAAAGCTGGCAAC 

+ 

208-228 6 

BCoV 

1503 

TAGATGAGATATCTGAACGACTGTCACC 

- 

112-119'’ 

BCoV 

1504 

CAGTCGTTCGCGAATTTCCGCAAAGAGTGGC 

+ 

221-240 6 

MHV-S 

1505 

CGGAAATTCGCGAACGACTGTCACCAAATAG 

- 

121-140 6 

MHV-S 

1509 

GACTGTAATCATGTTGTTAACACCAACCC 

+ 

127-155 6 

BCoV 

1510 

TCCGGAATTACACAGGGCAGGATTAAGGTC 

- 

175-198 6 

BCoV 

1555 

TGTACCTATGTGGAAAATAAC 

+ 

145-165 6 

MHV-S 

1556 

TTTTCCGGAATCACACAATTTCGGGTCAAG 

- 

193-210 6 

MHV-S 

1860 

CCACCATGTTTTTGCTTCTTAGATTTG 

+ 

1-22 6 

BCoV 

1861 

GGTCTAAGCATCATGCAGCC 

- 

1259-1278 6 

BCoV 

2060 

GAATTCCCACCATGGCGCGCACGGATGCA 

+ 

1-18 6 

DVIM 

2061 

TTATATTTTATGAAGGAGTTAATTT 

+ 

302-326'’ 

DVIM 

2063 

GATATCTTAAGCCTCATGCAATCTAAC 

- 

1276-1296'’ 

DVIM 

2078 

CTCCTTCATAAAATATAATTTGGGACCCTTCACCAGT 

- 

283-319 6 

DVIM 

2184 

TTTGGATCCACCATGTATTTGTTGCCGAAATTT 

+ 

1-21 6 

ECoV 

2185 

TTTCTGCAGCTAAGCTTCATGAAGCCTAGT 

- 

1252-1272'’ 

ECoV 


a Position with respect to the EToV Berne or MHV-A59 genome. 
b Positions as counted from the initiation codon of the respective nidovirus HE genes. 


Compton and Moore (AE091734). We therefore characterized the MHV- 
DVIM variant of the Compton laboratory (kindly provided by Dr. S. R. 
Compton). Reverse transcription-PCR products, together spanning a 
10.8-kb region comprising the 3'-end of ORFlb through the N gene, 
were cloned in pGEM-T Easy (Promega) and subjected to sequence 
analysis using the ABI Prism BigDye Terminators version 3.0 cycle 
sequencing kit (Applied Biosystems) and an ABI Prism 3100 genetic 
analyzer. The nucleotide sequence, determined in both orientations and 
for each region on at least two independent clones, was deposited in the 
EMBL data base (AY771998). Throughout, we use “DVIM” to refer to 
the Compton variant. 

Amino acid alignments were generated with Multalin, employing the 
Dayhoff symbol comparison table, and refined manually. Gaps were 
excluded. Phylogenetic trees were created by the neighbor-joining 
method, implemented in the program MEGA version 2.1 (44); the reli¬ 
ability of each interior branch was tested by bootstrap analysis with 
1000 resamplings. 

The number of synonymous substitutions per synonymous site (K s ) 
was estimated by sliding window analysis for overlapping 500-nucleo¬ 
tide gene segments with a 250-nucleotide step size with the program 
K-Estimator (45, 46), and D ss profiles were generated as described (20). 

Construction of HE Expression Vectors —Expression plasmid pBS- 
BRV-HE and reverse transcription-PCR cloning of the HE genes from 
torovirus strains B145, B150, P-MAR, P10, and P4 were described previ¬ 
ously (15, 20). The HE gene of MHV-S was reverse transcription-PCR- 
amplified from total intracellular RNA, isolated from infected cells. 2 
cDNA clones of the HE genes of MHV-DVIM, ECoV-NC99, and BCoV 
Quebec (47) were kindly provided by Drs. S. R. Compton, D. Brian, and 
A. A. Potter, respectively. The HE protein of BCoV Quebec differs from 
that of strain Mebus by a single substitution, Val 103 —> Leu, within a motif 
FYEGVN (residues 99—104), which is strictly conserved among nidovirus 
and influenza C virus O-acetylesterases; since this substitution might 
affect enzymatic activity, Leu 103 was changed back to Val via site-directed 
mutagenesis. All HE genes were inserted into pBluescript vectors (Strat- 
agene) downstream of the bacteriophage T7 RNA promoter, either as PCR 
amplicons generated with primer sets listed in Table I (ECoV, MHV- 
DVIM, and MHV-S) or by recloning restriction fragments from 
pcDNA3HE (BCoV) (47) or pGEM-T Easy vectors (all torovirus HEs). All 
constructs were sequenced. Inadvertent nucleotide substitutions were 
corrected in accordance with published sequences. 

Preparation of Sialic Acid Derivatives for Substrate Specificity As¬ 
says —Purified free Sia enriched for Neu5,9Ac 2 was obtained commer¬ 
cially (Sigma) and consisted of —17% Neu5Ac, 10% 5-iV-glycolylneura- 
minic acid, 2% 5-lV-acetyl-7-0-acetylneuraminic acid (Neu5,7Ac 2 ), 5% 
5-N-acetyl-8-0-acetylneuraminic acid (Neu5,8Ac 2 ), 34% Neu5,9Ac 2 , 
and 32% of a non-Sia contaminant. 


Sia enriched for Neu4,5Ac 2 was purified from horse serum glycopro¬ 
teins by acid hydrolysis (48). The mixture contained 3 mg of Sia/ml 
consisting of 69% Neu5Ac, 11% 5-iV-glycolylneuraminic acid, and 20% 
Neu4,5Ac 2 as determined by gas chromatography-mass spectrometry 
(see below). 5-iV-Acetyl-4,9-di-0-acetylneuraminic acid a-methyl glyco¬ 
side (uNeu4,5,9Ac 3 2Me) was synthesized by partial O-acetylation of 
Neu5Ac n-methyl glycoside (Sigma) as described previously (49). 

Heterologous Expression of Nidovirus HE Proteins and Substrate 
Specificity Assays —OST7-1 cells (—1 X 10 6 ) were infected with recom¬ 
binant vaccinia virus vTF7-3 and, at 1 h postinfection, transfected with 
2.5 pg of plasmid DNA as described (50). At 8 h postinfection, the cells 
were washed with phosphate-buffered saline and lysed in 200 pi of 
phosphate-buffered saline, 1% Triton X-100. Nuclei and cell debris were 
pelleted by centrifugation for 2 min at 14,000 rpm and 4 °C. Superna¬ 
tants were tested for sialate-O-acetylesterase activity with p-nitrophe- 
nyl acetate (pNPA) (Sigma) as a substrate (42). Lysates from mock- 
transfected vTF7-3-infected cells were included as negative controls. 
The specific enzymatic activity per ml of lysate, corrected for endoge¬ 
nous acetylesterase activity present in mock-transfected cells, was cal¬ 
culated as described (7). 

Substrate specificity assays were initially performed with equal 
amounts of specific enzymatic activity based upon hydrolysis of pNPA 
(7). However, as the activity of the acetylesterases toward pNPA proved 
of limited predictive value with respect to activity toward Sia deriva¬ 
tives, all assays were instead performed with fixed amounts of 
cell lysates. 

Substrate specificity of nidovirus HE proteins was determined in 100 pi 
of reaction mixtures containing 40 pi of cell lysate and 0.05 mg/ml 
Neu5,9Ac 2 , 0.3 mg/ml Neu4,5Ac 2 , or 0.1 mg/ml <xNeu4,5,9Ac 3 2Me in phos¬ 
phate-buffered saline, 0.5% bovine serum albumin. Incubation was for 1 h 
at 37 °C. Reactions were stopped by snap freezing and lyophilization. 

For assays with Neu5,9Ac 2 or Neu4,5Ac 2 substrates, Sias were deri- 
vatized with l,2-diamino-4,5-methylenedioxybenzene for 2.5 h at 50 °C 
(51) and analyzed by fluorimetric HPLC on a Cosmosil 5C18 ARII 
column (Nacalai Tesque; particle size 5 pm; 4.6 X 250 mm) with 
isocratic elution using water/methanol/acetonitrile (84:7:9, v/v/v; flow 
rate, 1 ml/min) and detection at excitation and emission wavelengths of 
373 and 448 nm, respectively. 

Assays with glycosidically bound Sias were performed essentially as 
described (42). Peak identities of the chromatograms were confirmed by 
MS analysis. Spectra were obtained on a Micromass Q-Tof Ultima 
Global mass spectrometer (Waters) with electrospray ionization and 
direct infusion. Dried samples were dissolved in 50% aqueous acetoni¬ 
trile containing 1% NH 4 0H and infused at a rate of 10 pl/min. Spectra 
were recorded in the negative ion mode for 5 min, using a scan time of 
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2 s and 0.1-s interscan delay, combined, and centered. Only singly 
charged ion species were observed. 

Enzymatic de-O-acetylation of aNeu4,5,9Ac 3 2Me was analyzed by 
gas chromatography-mass spectrometry (52). Lyophilizates were dis¬ 
solved in 40 pi of pyridine/hexamethyldisilazane/trimethylchlorosilane 
(5:1:1, v/v/v) and incubated for 45 min at room temperature to allow for 
trimethylsilylation of Sias. Samples (1 pi) were analyzed with a Fisons 
Instruments GC 8060/MD800 system (Interscience, Breda, The Neth¬ 
erlands) and an AT-1 column (30 m X 0.25 mm, Alltech, Breda, The 
Netherlands); the temperature program was 220 °C for 25 min, 6 °C/ 
min to 300 °C for 10 min; the injector temperature was 220 °C; detec¬ 
tion was done by electron impact MS with a mass range of 150-800 
mtz (53). 

Metabolic Labeling and Radioimmunoprecipitation Assay —Meta¬ 
bolic labeling of cells and subsequent radioimmunoprecipitation assay 
with rabbit-a-BToV BRV serum or with serum from a piglet naturally 
infected with PToV strain P-MAR (54) were performed as described 
(15, 55). 

RESULTS AND DISCUSSION 

Phylogenetic Relationships among Orthomyxo- and Nidovi¬ 
rus Sialate-O-acetylesterases —The HE proteins of toro- and 
coronaviruses and the influenza C virus HEF1 are evolutionary 
equidistant, sharing —30% sequence identity. Coronavirus 
HEs form disulfide-bonded homodimers (56); the oligomeric 
state of the torovirus HEs is unknown. Segments correspond¬ 
ing to HEF1 domains E and R and to what might be vestigial 
remnants of domain F (9) can be distinguished in the nidovirus 
HEs. The nidoviral HEs are membrane-anchored via a C-ter- 
minal hydrophobic domain, which in the primary sequence of 
HEF would correspond to the fusion peptide at the N terminus 
of HEF2. 

Most variation between and among orthomyxo- and nidoviral 
HEs occurs in domain R, and there is no obvious conservation 
of the HEF receptor-binding site in the nidovirus proteins. In 
domain E, the catalytic Ser-His-Asp triad (Ser 57 , His 355 , 
Asp 352 ) is strictly conserved. The HEF residues Gly 85 and 
Asn 117 , of which the NH group and side chain, respectively, 
together with the NH group of Ser 57 form the oxyanion hole (9), 
are maintained in most nidoviral HEs, but in those of the 
BToVs, Gly 85 is substituted by Ser (see Supplemental Fig. 1). 

In phylogenetic trees, the HE proteins of corona-, toro-, and 
influenza C viruses divide into three separate monophyletic clus¬ 
ters (Fig. 1; see Table II for abbreviations of virus names). This 
tree topology is consistent with the notion that the ancestors of 
torovirus, coronavirus, and influenza C virus acquired their HEs 
through independent heterologous recombination events (15) and 
that in each case the introduction of an HE gene into the viral 
genome occurred only once. Consequently, the coronavirus HE 
genes must all have evolved from one ancestral gene, which was 
acquired by a group II coronavirus predecessor after the split-off 
from HCoV-SARS. Similarly, the torovirus HEs all stem from a 
single ancestral O-acetylesterase. 

Whereas sequence variation among influenza C virus HEFs 
is limited, the corona- and torovirus HEs have diverged con¬ 
siderably (up to 40%). Four distinct torovirus HE lineages can 
be distinguished, represented by the HEs of BToV-BRV, BToV- 
B150, PToV-P4, and PToV P-MAR. A similar branching into 
four groups is seen for the coronavirus HEs (Fig. 1) (20). Pre¬ 
viously, only two lineages were recognized based upon phylo¬ 
genetic relationships and substrate specificity, with those of 
lineage A (represented by the HE of MHV-S and RCoV-SDAV) 
exhibiting sialate-4-O-acetylesterase activity and those of lin¬ 
eage B (represented by BCoV and HCoV-OC43) displaying 
sialate-9-O-acetylesterase activity (38). Our current analysis 
shows that, in fact, two additional lineages, C and D, can be 
distinguished (represented by the HEs of MHV-DVIM and 
MHV-2 and that of ECoV strain NC99, respectively), for which 
the substrate specificity remains to be determined (see below). 



C/MIY/4/93 

C/JHG/66 


Fig. 1. Phylogenetic relationships among orthomyxo- and 
nidovirus HE proteins. An alignment of torovirus, coronavirus and 
influenza C virus HE sequences (see Supplemental Fig. 1) was used to 
generate a neighbor-joining tree. Confidence values calculated by boot¬ 
strapping (1000 replicates) are indicated at the major branching points. 
Branch lengths are drawn to scale; the scale bar represents 0.1 amino 
acid substitution/site. The four different coronavirus lineages are indi¬ 
cated ( A-D). 

Table II 

Virus names and abbreviations 


Virus 

Strain/variant 

Abbreviation 

Coronavirus 

Bovine coronavirus 

Mebus 

BCoV 

Equine coronavirus 

NC99 

ECoV 

Porcine hemagglutinating 


PHEV 

encephalomyelitis virus 

Human coronavirus 

OC43 

HCoV-OC43 

Rat sialodacryoadenitis 


RCoV-SDAV 

coronavirus 

Puffinosis coronavirus 


PuCoV 

Mouse hepatitis virus 

S 

MHV-S 


A59 

MHV-A59 


JHM 

MHV-JHM 


DVIM 

MHV-DVIM 


2 

MHV-2 

Torovirus 

Bovine torovirus 

Breda 

BToV-BRV 


B145 

BToV-B145 


B150 

BToV-B150 

Porcine torovirus 

Markelo 

PToV-P-MAR 


P10 

PToV-PIO 


P4 

PToV-P4 

Orthomyxovirus 

Influenza C virus 

Miyagi/4/93 

C/MIY/4/93 


Johannesburg/66 

C/JHG/66 


Torovirus HEs: PToV and BToV HE Proteins Differ in Their 
Preference for O-Acetylated Sialic Acid Subsets —To gain more 
insight into the evolution of the nidoviral O-acetylesterases, we 
studied the substrate specificities of the HEs of coronaviruses 
MHV-DVIM and ECoV and of representatives of each of the 
four torovirus lineages. Since toroviruses, with Berne virus as 
the sole exception, cannot be propagated in tissue culture, the 
various HE proteins were produced in the vaccinia virus-based 
vTF7-3 expression system. Successful expression of each viral 
acetylesterase was confirmed by radioimmunoprecipitation as- 
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Fig. 2. Heterologous expression of nidovirus HE proteins. A, 

detection of vTF7—3-expressed nidovirus HE proteins by radioimmuno- 
precipitation. vTF7-3-infected cells, (mock-)transfected with HE ex¬ 
pression plasmids were metabolically labeled. Lysates were subjected to 
a radioimmunoprecipitation assay with a rabbit antiserum raised 
against BToV BRV ( aBRV) or with serum from a piglet, naturally 
infected with PToV P-MAR (aP-MAR). Samples were analyzed in SDS- 
polyacrylamide gels. Molecular size markers are given in kDa. B, de¬ 
tection of vTF7—3-expressed nidovirus HE proteins by acetylesterase 
assay. Lysates from vTF7—3-infected, (mock-)transfected cells were har¬ 
vested at 8 h postinfection. Samples (10 pi) were tested for acetylester¬ 
ase activity by incubation with 1 mM pNPA in phosphate-buffered 
saline in 1-ml reaction volumes. Cumulative specific hydrolysis of 
pNPA (in nmol), measured at 405 nm and corrected for background 
acetylesterase activity in mock-transfected cells (y axis), is plotted 
against reaction time (x axis). 


say (BToV strain BRV and B145 and PToV strain P-MAR) 
and/or by demonstrating enzymatic activity in cell lysates with 
pNPA as a substrate (Fig. 2 and Table III). Note that antisera 
against BToV-BRV also detected the HE of BToV-B145 but not 
that of BToV-B150 (Fig. 2A), illustrating that the two BToV HE 
lineages (Fig. 1) are antigenically distinct. 

Substrate preferences of the nidoviral HEs were studied with 
mixtures of purified free mono-O-acetylated Sias, enriched for 
Neu4,5Ac 2 or Neu5,9Ac 2 . In agreement with published obser¬ 
vations (22, 41, 42), vTF7-3-expressed HEs of MHV-S and 
BCoV-Mebus, taken along as model enzymes, preferentially 
removed 4-0- and 9-O-acetyl groups, respectively (Fig. 3A and 
Table III). The HEs of the porcine torovirus strains P-MAR, 
P10, and P4 all displayed a clear preference for Neu5,9Ac 2 (Fig. 
3 B and Table III). Bovine torovirus HEs displayed only mod¬ 
erate enzymatic activity toward either Sia substrate; although 
Neu5,9Ac 2 was hydrolyzed preferentially, Neu4,5Ac 2 was 
cleaved as well (Fig. 3C and Table III). 

The fact that free Sias mainly adopt the j3- (>93%) instead of 
the a-configuration, present in glycosidically bound Sias, posed 
a potential caveat, since substrate anomeric configuration may 
well influence enzymatic cleavage (11). Therefore, we also as¬ 
sessed enzyme specificity with aNeu4,5,9Ac 3 2Me, a synthetic 
substrate with an a-aglycon group. The HEs of MHV-S and 
BCoV-Mebus again preferentially cleaved 4-0- and 9-O-acetyl 


Table III 

Hydrolysis of pNPA and of mono- and di-O-acetylated Sias 


Virus 

pNPA 

Mono-0-Sia a 

Di-O-Sia 6 

4-0 

9-0 

4-0 

9-0 


milliunits/ml 

% 

% 

% 

% 

MHV-S 

30 

96 

i 

100 

16 

BCoV 

119 

6 

56 

0 

100 

MHV DVIM 

50 

2 

50 

0 

95 

ECoV 

120 

4 

58 

0 

100 

BToV BRV 

39 

13 

14 

0 

28 

B145 

43 

14 

28 

0 

65 

B150 

16 

7 

25 

0 

98 

PToV P-MAR 

31 

6 

96 

0 

100 

P10 

23 

4 

84 

0 

100 

P4 

7 

0 

23 

0 

29 


“ Purified free mono-O-acetylated sialic acids Neu4,5Ac 2 (4-0) and 
Neu5,9Ac 2 (9-0). 
b ciNeu4,5,9Ac 3 2Me. 


groups, respectively. When tested with aNeu4,5,9Ac 3 2Me, the 
HEs of both PToV and BToV exclusively hydrolyzed the 9-O- 
acetyl group (Fig. 4 and Table III). 

These experiments tentatively identified the toroviral HEs 
as sialate-9-O-acetylesterases. We noted, however, a consistent 
difference between BToV and PToV enzymes; although the 
HEs of BToV strains BRV and B145 and those of PToV strains 
P-MAR and P10 hydrolyzed the synthetic substrate pNPA to 
similar extents, the BToV enzymes were reproducibly less ac¬ 
tive toward free Sia substrates (Table III). This raised the 
question of whether a difference in substrate specificity be¬ 
tween the PToV and BToV O-acetylesterases might yet exist. 
We therefore compared the reactivity of these enzymes toward 
bovine submaxillary mucin (BSM), a model substrate, which 
contains a mixture of glycosidically bound Sias (Table IV). The 
HE of PToV P-MAR hydrolyzed 9-0-acetyl-5-A r -glycolylneura¬ 
minic acid, Neu5,9Ac 2 , Neu5,8Ac 2 , and possibly also Neu5,7Ac 2 
but showed little reactivity toward 5-N-acetyl-7(8),9-di-0- 
acetylneuraminic acid (Neu5,7(8),9Ac 3 ) (Fig. 5A and data not 
shown). In contrast, the HE of BToV-BRV exhibited only mod¬ 
est activity toward mono-O-acetylated Sias but converted 
Neu5,7(8),9Ac 3 very efficiently (Fig. 5A), apparently to 
Neu5,7Ac 2 and Neu5Ac (Fig. 5 B). Similar observations were 
made for the HEs of BToV strains B145 (not shown) and B150 
(Fig. 5A). None of the enzymes showed any reactivity when 
tested with horse serum glycoproteins, a model substrate con¬ 
taining glycosidically bound Neu4,5Ac 2 (not shown). Our find¬ 
ings suggest that the PToV and BToV HEs indeed differ in 
substrate specificity, displaying distinct preferences for Sia 
subsets with either a mono -O- or di-O-acetylated glycerol side 
chain, respectively. 

Coronavirus HEs: The HE Protein of ECoV-NC99 Is a 
Sialate-9-O-acetylesterase —The HE proteins of MHV strains 
DVIM/-2 (lineage C) and ECoV-NC99 (lineage D) are distinct 
from each other and from those of the coronavirus HE lineages A 
and B (Fig. 1). The HE of ECoV is most closely related to the 
sialate-9-O-acetylesterases of the BCoV/HCoV-OC43 cluster, but 
since 4-0-acetylated Sias are prominently expressed in the horse 
(2), we entertained the possibility of ECoV HE being a sialate-4- 
O-acetylesterase. However, in cleavage assays with the free Sia 
substrates Neu4,5Ac 2 , Neu5,9Ac 2 , and «Neu4,o,9Ac 3 2Me, ECoV 
HE was indistinguishable from BCoV HE and preferentially 
cleaved at C9 (Figs. 3 and 4; Table III). Also, when tested with 
glycosidically bound Sias of BSM, ECoV HE preferentially de-O- 
acetylated Neu5,9Ac 2 (not shown). These observations identify 
ECoV-NC99 HE as a classical sialate-9-O-acetylesterase. 

Coronavirus HEs: The HE Protein MHV-DVIM Is Not a 
Sialate-4-O-acetylesterase —We anticipated MHV-DVIM, like 
all other rodent coronaviruses studied so far, to encode a 
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Fig. 3. Reactivity of heterologously expressed nidovirus HE 
proteins toward purified free mono-O-acetylated Sias. Free Sia 
preparations, enriched for Neu5,9Ac 2 ( peak 5, upper panels ) or 
Neu4,5Ac 2 ( peak 6, lower panels), but also containing 5-iV-glycolyl- 
neuraminic acid ( peak 1 ), Neu5Ac (peak 2), Neu5,7Ac 2 (peak 3), and 
Neu5,8Ac 2 (peak 4) and a non-Sia compound (peak 7) were incubated 
with vTF7—3-expressed HE proteins of coronaviruses (A), porcine toro- 
viruses (B), and bovine toroviruses (C). Sias were derivatized with 
l,2-diamino-4,5-methylenedioxybenzene and analyzed by fluorimetric 
reversed-phase C 18 HPLC. In the chromatograms, fluorescence at 448 
nm (y axis) is plotted against retention time (x axis). 


sialate-4-O-acetylesterase. However, in contrast to MHV-S HE, 
DVIM HE did not cleave the 4-0-acetyl groups from the free 
Sia substrates Neu4,5Ac 2 and nNeu4,5,9Ac 3 2Me (Figs. 3 and 4) 
nor from glycosidically bound Neu4,5Ac 2 present in horse se¬ 
rum glycoproteins (not shown). Instead, DVIM HE hydrolyzed 
the 9-0-acetyl moieties from free Neu5,9Ac 2 and from 
nNeu4,5,9Ac 3 2Me. Puzzlingly, the enzyme showed little reac¬ 
tivity toward any of the glycosidically bound Sias of BSM (Fig. 
5). Substrate specificity of DVIM HE might be restricted to a 
subset of 9-O-acetylated Sias not present in BSM (e.g. efficient 
cleavage might require additional Sia modifications or might 
be linkage-dependent). Recently, we have found that treatment 


Mock MHV-S BCoV MHV-DVIM ECoV 
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Fig. 4. Reactivity of heterologously expressed nidovirus HE 
proteins toward aNeu4,5,9Ac 3 2Me. The synthetic Sia derivative was 
incubated with vTF7-3-expressed nidoviral HE proteins, derivatized by 
trimethylsilylation, and analyzed by gas chromatography-mass spec¬ 
trometry. The resulting total ion current chromatograms are presented. 
The following Sias were identified by mass spectrometry: uNeu5Ac2Me 
(peak 1), uNeu4,5,9Ac 3 2Me (peak 2), aNeu5,9Ac 2 2Me (peak 3), 
aNeu4,5Ac 2 2Me (peak 4)\ peaks 2' and 5 represent non-Sia compounds. 



Table IV 

Sialic acid content of BSM 

Sialic acid Percentage of total 

sialic acid 


5-V-glycolylneuraminic acid 

12 

Neu5Ac 

22 

Neu5,7Ac 2 

13 

9-0-acetyl-5-V-glycolylneuraminic acid 

9 

Neu5,8Ac 2 

5.5 

Neu5,9Ac 2 

28 

Neu5,7(8),9Ac 3 

10.5 


of chicken erythrocytes with DVIM HE destroyed BCoV recep¬ 
tors (i.e. Neu5,9Ac 2 ), thus preventing BCoV-mediated hemag¬ 
glutination; conversely, treatment of rat erythrocytes with 
BCoV HE abrogated agglutination by MHV-DVIM. 4 The com¬ 
bined observations strongly suggest that MHV-DVIM binds to 
(a subset of) 9-O-acetylated Sia receptors and that its HE is a 
sialate-9-O-acetylesterase. We conclude that murine coronavi¬ 
ruses occur in the field in two types, which differ with respect 
to Sia receptor usage and HE substrate specificity. 

Evidence for Horizontal Transfer of HE Genes among Coro¬ 
naviruses —It is important to stress that the phylogenetic rela¬ 
tionships in Fig. 1 reflect the evolutionary history of the HE 
modules and not necessarily that of the viruses from which 
these genes were derived, since nidoviruses frequently engage 
in homologous RNA recombination (57). Indeed, there is evi¬ 
dence for recurrent horizontal transfer of HE genes among 
toroviruses; ancestors of bovine and porcine field strains ap¬ 
parently replaced their original HE sequences via intertypic 
RNA recombination with hitherto unidentified toroviruses, 


4 A. L. W. van Vliet, M. Langereis, and R. J. de Groot, unpublished 
observations. 


Downloaded from http://www.jbc.org/ at University of California, Santa Cruz on March 11, 2015 

























































Evolution and Substrate Specificities of Nidovirus HE Proteins 


6939 


A 

-O- mock -A- P-MAR HD- BRV -o- B150 DVIM 




Fig. 5. Reactivity of heterologously expressed nidovirus HE 
proteins toward glycosidically bound 9-O-acetylated sialic acids 

of BSM. A, BSM was incubated with vTF7—3-expressed nidoviral HE 
proteins for 0, 1, 3, and 18 h at 37 °C. After acid hydrolysis and 
derivatization with l,2-diamino-4,5-methylenedioxybenzene, Sias were 
analyzed by fluorimetric reversed-phase HPLC. Graphs depict the en¬ 
zymatic hydrolysis of Neu5,7(8),9Ac a (upper panel) and Neu5.xAc 2 (i.e. 
Neu5,7Ac 2 , Neu5,8Ac 2 , and Neu5,9Ac 2 combined) (lower panel) over 
time. B, sialic acid composition of BSM after (mock) treatment for 0, 1, 
and 3 h with BToV-BRV HE. In the bar graphs , relative amounts of 
Neu5,7(8),9Ac 3 (5,7 (8),9), Neu5Ac (5), Neu5,7Ac 2 (5,7), Neu5,8Ac 2 (5,8), 
and Neu5,9Ac 2 (5,9) are depicted as percentages of total BSM 
Sia content. 

thus giving rise to the BToV B150 and (presumably) the PToV 
P-MA.R/P10 lineages (20). The fact that in the coronavirus 
branch, the HEs of the MHV strains DVIM and MHV-2 cluster 
separately from those of other MHV strains (JHM, S, and A59) 
and of RCoV-SDAV (Fig. 1) prompted us to consider the possi¬ 
bility that, also during coronavirus evolution, an exchange of 
HE genes had occurred. 

As the sequence data available for MHV-DVIM are incom¬ 
plete and conflicting (also see “Materials and Methods”), we 
determined the nucleotide sequence for the 3' one-third of the 
genome for the DVIM stock obtained from the Compton labo¬ 
ratory. Phylogenetic trees generated per gene revealed incon¬ 
gruencies, indicative for recombination (Fig. 6A). In the HE 
gene, MHV strains DVIM and MHV-2 grouped separately from 
the other rodent coronavirus strains (A59, JHM, and RCoV- 
SDAV). In the S gene, MHV-DVIM and MHV-2 were also 


distinct from strains A59 and JHM but now grouped together 
with RCoV-SDAV. Finally, in the surrounding genes (ORFlb, 
ns2, E, M, and N), all rodent coronaviruses were closely related 
and in essence formed one single cluster. 

Further support for recombination was obtained by analysis 
of the distribution of synonymous substitutions (D ss ; Fig. 6 B). 
D ss profiles based upon pairwise comparison of MHV-DVIM 
and -2 with JHM and A59 showed that within the HE and S 
genes, K s values (i.e. the estimated number of synonymous 
substitutions per synonymous site) (46) are >6-fold higher 
than in the flanking genes. Comparison of RCoV-SDAV versus 
MHV-A59 yielded roughly constant K s values across the ns2, 
HE, E, M, and N genes, but in the S gene, the number of 
synonymous substitutions per synonymous site was 5-fold 
higher. Most interestingly, pairwise comparison of MHV-DVIM 
and MHV-2 with RCoV-SDAV revealed that these three viruses 
are very closely related. RCoV-SDAV only differs from these 
MHV strains in the HE gene as indicated by a 5-fold increase in 
the K s values. Surprisingly, however, MHV-DVIM and MHV-2 
also seem to be recombinant with respect to each other as 
indicated by an abrupt 3-fold rise in K s values in the HE gene. 

The combined observations indicate that the divergence of 
the rodent coronaviruses was indeed accompanied by multiple 
recombination events during which both HE and S sequences 
were replaced. Apparently, related strains as well as hitherto 
unidentified coronaviruses served as donors. Our current data 
do not allow conclusions as to which of the MHV strains rep¬ 
resent parental types and which represent recombinant prog¬ 
eny. The exchange among rodent coronaviruses of S and HE 
genes is highly reminiscent of the reassortment of HA and NA 
genome segments, which gives rise to antigenic shifts in influ¬ 
enza A viruses. We propose that, in analogy, the observed 
murine coronavirus diversity, brought about by RNA recom¬ 
bination, is a direct consequence of immune pressure for anti¬ 
genic variation. In the case of the MHV HE proteins, the genetic 
exchanges would, however, not only have resulted in an antigenic 
shift, but, at least on one occasion, also in a shift in substrate 
preference from 4-0- to 9-O-acylated Sias or vice versa. 

The Evolution of Nidoviral Acetylesterases: A Complex Tale 
of Divergence and Recombination —In comparison with other 
nidovirus structural proteins, the HEs have received little at¬ 
tention, probably because they are dispensable for replication 
in vitro and, in the case of the coronaviruses, do not belong to 
the standard protein repertoire. Yet, the mere fact that in the 
course of evolution, two groups of nidoviruses independently 
acquired an HE gene indicates that these virion-associated 
enzymes must provide a huge selective advantage during the 
natural infection (15). Our current observations strongly sup¬ 
port this view: under field conditions, HE genes may be ex¬ 
changed via homologous RNA recombination both among coro¬ 
naviruses and among toroviruses, but HE expression per se is 
strictly maintained. This is in stark contrast with the loss of 
HE expression, which occurs during virus propagation in tissue 
culture cells (13, 15, 26). 2 

The evolution of the coronavirus and torovirus HE proteins has 
yielded a variety of sialate-O-acetylesterases. Current diversity 
seems to have arisen not only from antigenic drift but also as a 
result from adaptations to new hosts and/or new target tissues. 
Under the assumption that HE substrate preference reflects Sia 
receptor usage, nidoviruses apparently employ Neu5,9Ac 2 , 
Neu4,5Ac 2 , and Neu5,7(8),9Ac 3 as initial attachment factors. As 
first reported here, HEs with a preference for the latter Sia 
subtype occur in bovine toroviruses, both in the BRV/B145 and 
B150 lineages. The latter observation is of particular interest, 
since the B150-type BToVs are recombinant viruses, which ap¬ 
parently originated from an exchange of HE sequences between 
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Fig. 6. Evidence for multiple exchanges of S and HE sequences 
during the radiation of rodent coronaviruses. A, rooted neighbor¬ 
joining trees depicting the phylogenetic relationships among rodent 
coronavirus strains. Trees were constructed for the 3'-most 1476 resi¬ 
dues of ORFlb and for the ns2, HE, S, E, M, and N genes with BCoV 
serving as an outgroup. Confidence values are indicated at the major 
branching points. Branch lengths are drawn to scale; the scale bar 
represents 0.05 nucleotide substitution per site. B, distribution of syn¬ 
onymous substitutions ( D ss ) in the lb, ns2, HE, S, E, M, and N genes of 
rodent coronavirus variants. Sequences were compared pairwise by 
sliding window analysis; the number of synonymous substitutions per 
synonymous site (K s ) was estimated for overlapping 500-nucleotide 
gene segments with a 250-nucleotide step size. D SB profiles were gene- 


a B 145-type ancestor and an as yet unidentified torovirus donor 
(20). This raises the question of whether the donor virus already 
encoded a sialate-7 (8)9-di-0-acetylesterase or if the substrate 
preference of the newly acquired HE was adjusted after the 
horizontal gene transfer. Porcine toroviruses at least display a 
preference for Neu5,9Ac 2 receptors. 

The divergence of the coronavirus HEs produced four distinct 
lineages. Whereas the HEs of lineage A preferentially hydro¬ 
lyze 4-O-acetylated Sias, those from the three other lineages 
de-O-acetylate Sias preferentially at C9. Based upon phyloge¬ 
netic tree topologies (Fig. 1 and data not shown), we propose 
that the ancestral coronavirus HE was a sialate-9-O-acetyles- 
terase, which was converted into a sialate-4-O-acetylesterase 
presumably through the accumulation of mutations in the E 
domain. The substrate binding site of the coronavirus sialate- 
4-O-acetylesterases must differ essentially from that of their 
9-O-specific homologues; bringing the 4-0- instead of the 9-0- 
acetyl group in close proximity of the principal nucleophile 
would require a rotation of the Sia with respect to the catalytic 
site. In which coronavirus and host the transition in substrate 
preference took place is not known; nor is the direction of the 
various recombinational exchanges that occurred during sub¬ 
sequent divergence of the rodent coronaviruses. Possibly, the 
sialate-4-O-acetylesterase arose in an as yet unknown corona¬ 
virus to be introduced into the genome of a DVIM-type MHV 
via horizontal gene transfer. Irrespective of the precise course 
of events, it has led to the natural occurrence of two types of 
closely related murine coronaviruses. These variants differ 
with respect to substrate preference of their HE proteins, im¬ 
plying that they use different Sia receptors as attachment 
factors to initiate their infections in the mouse. It is of note that 
both 4-0- and 9-O-acetylated Sias as well as less common Sias 
like 5-lV-acetyl-9-0-lactylneuraminic acid are present in mouse 
colon and trachea and in virtually all other target organs of 
MHV (58). 5 The consequences of Sia receptor specificity for host 
cell tropism and disease outcome may be considerable (59) and 
remain to be explored. 

The HE-expressing coronaviruses seem to have adopted a 
lifestyle that in many ways resembles that of influenza A 
viruses. There are parallels with respect to virion composition 
and host cell attachment in that Sia receptor binding, and RDE 
activities are performed by two different surface glycoproteins 
(HA and S versus NA and HE, respectively). These envelope 
proteins occur in multiple subtypes, which differ in antigenic 
properties and in Sia receptor specificity casu quo substrate 
preference. Antigenic shifts may occur by exchanging the genes 
(or parts thereof) for the envelope glycoproteins, in concert or 
individually, through reassortment or RNA recombination. 
However, ultimate success of the recombinant offspring would 
depend not only on its potential to escape host and population 
immunity but also on proper matching of receptor-binding and 
RDE activities. Depending on their host species, influenza A 
virus strains differ in their preference for particular Sia sub¬ 
sets (e.g. avian strains preferentially bind to sialoglycoconju- 
gates containing «2,3-linked Sias, whereas human strains use 
a2,6-linked Sia-receptors). A mismatch between receptor pref- 


5 A. Rinninger, C. Richet, A. Pons, G. Kohla, R. Schauer, H.-C. Bauer, 
R. Vlasak, J.-P. Zanetta, manuscript in preparation. 


rated by plotting calculated K s values against nucleotide position. As a 
reference, the genes for lb, ns2, HE, S, E, M, and N, drawn to scale, are 
depicted as boxes (top). C, simplified schematic representation of the 
recombinant genome composition of rodent coronaviruses. Genes are 
represented as boxes. Related S and HE genes are color-coded. The 
genome structure indicated by a question mark has not yet been iden¬ 
tified in the field. 
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erence and RDE substrate specificity brought about by reas¬ 
sortment causes a decrease in viral fitness, which can be re¬ 
stored by compensatory mutations in HA, NA, or both proteins 
(5). A similar functional link might exist between the nidovirus 
S and HE proteins. It is of interest that in rodent coronavi- 
ruses, two major types of S proteins can be distinguished 
(JHM/-S versus DVIM/RCoV; Fig. 6, A and B ), that (parts of) S 
and HE genes are exchanged individually, and that three of 
four possible S-HE combinations have been found so far (Fig. 
60. Intriguingly, the rodent coronavirus strains MHV-DVIM 
and RCoV-SDAV express S proteins, which are closely related 
(90% identity) but have HE proteins with different substrate 
specificities. Future studies should (i) determine the Sia recep¬ 
tor preference of these viruses, (ii) determine whether also in 
these viruses the S proteins are the major Sia-binding proteins, 
and (iii) if so, establish whether the receptor-binding and RDE 
activities of their respective S and HE proteins match. 

Our findings provide a framework for further structure-func¬ 
tion analysis of nidovirus sialate-O-acetylesterases, which 
should aim at the elucidation of the molecular basis for sub¬ 
strate specificity and explore the potential of HE-directed an¬ 
tiviral drugs. Moreover, the HEs may serve as novel tools to 
study sialic acid biology. 
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